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^ We test numerical filamentation models against experimental data about 
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quantitative agreement without the need of adjustable parameters. ©2012 
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Laser filamentation [1-4] is a propagation regime typical of high-power lasers. It stems from 
a dynamic balance between self-focusing and self-defocusing non-linearities of different or- 
ders, i.e., different intensity dependences. Self-focusing by the Kerr effect is well established. 
Conversely, the relative contributions to defocusing are still discussed. The two main ef- 
fects nowadays considered in gases are the laser-generated plasma and the saturation of the 
medium polarisability under strong- field illumination [5-7], empirically described as a power 
series of the intensity and referred to as higher-order Kerr effect (HOKE). The latter have 
been first introduced as freely adjustable parameters [8,9], before experimental values were 
reported [10, 11]. 

Strong effort have been dedicated to reach a quantitative agreement between numerical 
models of non-linear pulse propagation and experimental data about filamentation [12-15]. 
In that purpose, a common approach consists in tweaking simulation inputs like the non- 
linear index or the ionization rates, and/or the experimental parameters such as the beam 
energy, pulse duration, beam diameter or profile, until a satisfactory match with the ex- 
periments is achieved. The remaining discrepancies are generally attributed, among others, 
to the difficulty of the experimental measurements inside the filaments, and the associated 
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Fig. 1. Experimental setup 



uncertainties. However, the introduction of new techniques, e.g. for measuring the peak in- 
tensity from the ratio of nitrogen fluorescence lines [16,17] enhances the precision and renew 
the challenge for theoretical models to precisely match those data. 

Here, we provide experimental electron density measurements over a wide range of pulse 
durations, and compare these experimental results with numerical simulations. We show that 
considering the HOKE in the model allows a good agreement between experimental results 
and simulations without any parameter adjustment. Such agreement is also confirmed when 
comparing the model including the HOKE [18] with previously published data about the 
filament peak intensity [17] and electron density [19]. 

The filament ionization was characterized as a function of the pulse duration by using 4 mJ 
laser pulses centered at 800 nm, with 100 Hz repetition rate and 3 cm initial beam diameter. 
The slightly diverging beam was focused by an f=2.8 m lens to form a 20 cm long filament 
in air. The pulse duration was varied from 80 fs to 1.6 ps by detuning the grating compressor 
of the chirped pulse amplification (CPA) chain. The charge is measured by propagating the 
filaments between a pair of parallel electrodes (lxl cm; 1 cm spacing) under 1 kV bias: 
the transient current through the circuit is proportional to the electron density generated by 
the filaments [20], although with an arbitrary calibration factor (See Figure 1). This charge 
measurement was double-checked by sono metric measurements [21]. 

To get information about the role of the HOKE, we performed numerical simulations con- 
sidering a linearly polarized incident electric field at a wavelength Ao=800 nm with cylindrical 
symmetry around the propagation axis z. According to the unidirectional propagation pulse 
equation (UPPE) [22], the scalar envelope e(r, t, z) (defined such that \e(r, z, t)\ 2 = I(r, z, t), 
I being the intensity) evolves in the frame traveling at the pulse velocity according to [23]: 



where c is the velocity of light in vacuum, oj is the angular frequency, k(cu)—n(uj)cj / c, k' its 
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Fig. 2. Effect of pulse duration on the electron density in filaments. Each numerical curve is 
normalized independently to optimally match the experimental results. 

derivative at ujq=2ixc/ \ Q) n(u) is the linear refractive index at the frequency u, k± is the 
spatial angular frequency. Pnl is the nonlinear polarization, J is the free-charge induced 
current and a is the nonlinear losses induced by photo-ionization. / denotes simultaneous 
temporal Fourier and spatial Hankel transforms of function /. The non-linear polarization is 
evaluated as Pnl = n 2j\^ + An Raman j e, where the n 2 j are the j th -order nonlinear re- 
fractive indices (n 2 =1.2x 10~ 7 cm 2 /TW, n 4 =-1.5xl0- 9 cm 4 /TW 2 , n 6 =2.1xl(r 10 cm 6 /TW 3 , 
72 8 =-8xl0~ 12 cm 8 /TW 4 [10]) and AnR aman is the Raman induced refractive index change 
evaluated by solving the rotational time-dependent Schrodinger equation of both N 2 and 
O2 in the weak field regime. Alternatively, the model is truncated to the third-order non- 
linearity (i.e. the n 2 |e| 2 £ term) to disregard the contribution of the HOKE. The current is 
evaluated as J = ^-(^c + 2w)/5£/(V 2 + uj2 )i where e and m e are the electron charge and mass 
respectively, e is the vacuum permittivity, u e is the effective electronic collisional frequency, 
and p is the electron density. Finally, a = J2i=o 2 n 2 Wi(\£\ 2 )UiPat t i/ (2\e\ 2 ) , p at ,j is the density 
of molecules of species i, Wi(\e\ 2 ) is their photoionization probability modeled by the PPT 
formulation, with ionization potential C/j. 

The propagation dynamics of the electric field is coupled with the electron density p, 
calculated as [2] 



where (3 are the electron recombination rate and <Xj are the inverse Bremsstrahlung cross- 
sections with species i, also accounting for avalanche ionization. Note that our model relies 
only on published values [4,10] and contains no adjustable parameter. For each experimental 
situation, we directly used the experimental parameters, without any tweaking. Unless oth- 
erwise specified, Gaussian pulse shapes have been considered, both spatially and temporally. 
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Fig. 3. Experimental [17] and simulated peak intensity in filaments generated by 42 fs pulses 
of 1 cm diameter, focused with (a) / = 1 cm and (b) /=50 cm. 

The measured electron density rose very slowly (less than a factor of 4) when the pulse 
duration was decreased from 1.6 ps to 80 fs, although the peak power drastically increased by 
a factor of 20 (Figure 2). This might seem unexpected when considering that the ionization 
rate scales with I 8 in oxygen [24,25]. This relative stability can however be understood by 
considering that the free electrons need several picoseconds to tens of ps to recombine, so 
that they accumulate more efficiently during longer pulses, an effect that partly balances the 
lower ionization rates. While both models reproduce this slow variation for pulses below a 
few hundreds of femtoseconds, the full model seems slightly more accurate for shorter pulses. 

To get a more direct comparison between the two models, we successively tested them 
against previously published quantitative experimental data about the peak intensity [17] 
and electron density [19] measured in filaments by different groups and techniques, so as to 
avoid individual artifacts or systematic errors to affect our conclusions. 

As shown in Figure 3, the peak intensity of filaments generated by focused 42 fs pulses of 
1 cm diameter at 1/e 2 (i.e, 5.9 mm FWHM) measured via the ratio of nitrogen fluorescence 
lines [17] is much better reproduced by the full model than by the truncated one. The latter, 
in which the HOKE defocusing terms are disregarded, overestimates the intensity by a factor 
of typically 2-3. 

The electron density as measured in [19] provides a much more sensitive variable, be- 
cause its high-order dependence on intensity magnifies the experiment dynamics, allowing 
more reliable comparisons. The intensity of the nitrogen fluorescence constitutes a classical 
measurement of the electron density [26]. Consistent with the differences in the predicted 
intensities, the full model always yields a 4- to 30-times lower electron density than the 
truncated one. The match with experimental data [19] is perfect for a parallel beam cor- 
responding to typical filamentation conditions (Figure 4a). However (Figure 4b-d), even 
the full model tends to overestimate the electron density produced by focused beams. This 
discrepancy could be due to the recently suggested overestimation [7] of the free electron 
density by the PPT model [4, 24] , or to optical aberrations and astigmatism in experiments, 
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Fig. 4. Experimental [19] and simulated electron density as a function of focusing and incident 
power, for a 45 fs pulse of initial beam diameter of 8.4 mm (full width at 1/e 2 ), for (a) a 
free-propagating (parallel) beam; (b) / = 380 cm; (c) / = 100 cm; (d) / = 50 cm. 

which are not included in the model and tend to decrease the peak intensity in the waist 
region [27] and reduce filament length and strength [28] . Besides, this increasing geometrical 
constraint explains the decrease of the difference between the two models for tighter focusing. 
Therefore, the most concluding conditions to test filamentation models correspond the less 
focused ones, i.e. the larger / number. 

Our results show that considering the saturation of the medium polarization in the strong- 
field regime (aka HOKE) improves the quantitative modeling of filamentation, especially for 
collimated beams for which the geometrical constraint does not interfere with the self-guiding 
process. Such update of the standard laser filamentation model is easy to implement, based on 
both experimental measurement [10] and quantum mechanical justifications [5-7]. It will in 
particular impact the determination of the optimal conditions for the potential atmospheric 
applications [29,30] of laser filamentation, like rainmaking [31], or lightning control [32]. 
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